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I—ON THE PROPAGATION OF A WAVE IN AN 
ELASTIC MEDIUM. 


Tue following investigation of the equation for the propagation of 
a wave in an elastic medium, may be considered as supplemental to 
the two articles on the Undulatory Theory in the first and second 
Numbers of this Journal: though perhaps it ought more properly 
to have preceded them. But as we do not pretend to offer a com- 
plete treatise on the subject, the order of the articles is not a matter 
of much consequence. 

1. Light, according to the wave theory, consists of undulations 
propagated in an elastic medium. From the laws of the interference 
of polarized light, it appears that the constituent vibrations are 
transverse to the direction of propagation ; that polarized light con- 
sists of vibrations in one direction, or perpendicular to one plane, 
which in Fresnel’s theory is the plane of polarization, and in 
Cauchy’s theory is perpendicular to the plane of polarization ; and 
thirdly, that there are no vibrations, or at least none capable of 
producing the sensation of light, in the direction of propagation. 
Fresnel assigns as the reason for this, that the resistance to com- 
pression of the medium is so great, that if a disturbance takes place 
motion is propagated to a great distance before the disturbed mole- 
cule returns to rest. Hence, a wave of great breadth will be pro- 
pagated, but of feeble intensity, as the wis viva remains constant. 


2. To this we may add, that the possibility of the propagation 


breadth can be propagated: if we examine our equations, we see 
that this depends upon the disturbance being a function of vt—a; 
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and the partial differential equation which gives this form of solu- 
tion is obtained, in the simplest case, namely that of motion in one 
direction, by supposing the elastic force proportional to the con- 
densation, and the condensations to be infinitesimal. Without 
these restrictions the equations cannot be solved; but from several 
considerations it seems probable, that a particle once disturbed does 
not immediately come to rest, and that the waves increase in 
breadth. We may see examples of this in the waves propagated 
when a stone is thrown into still water. But in water the velocity 
of a wave depends on its breadth; and thus, though the breadths of 
the waves increase, they do not interfere with each other. 

3. It appears probable, from the change in the character of 
sounds caused by distance, the loss of sharpness and distinctness, 
that the breadth of a wave of air increases. If we suppose the 
same to be the case with those waves of light which consist of vi- 
brations in the direction of propagation, but to a greater degree 
than in sound, we shall have, in the first place, a rapid diminution 
in the vis viva of each wave ; and also, as their velocities are nearly 
independent of their breadths, each wave will overtake its prede- 
cessor and be overtaken by its successor, and there will be an 
interference which may speedily obliterate all traces of the disturb- 
ance. 

1. We have now to consider that complex system of disturbances 
which constitutes an undulation. 

We shall suppose the molecules of the media which we consider 
to be so arranged, that the three rectangular directions, or the three 
axes of elasticity, have the same direction at every point. These 
we shall take for coordinate axes, and shall suppose that the elas- 
ticity called into play is greatest for a displacement parallel to the 
axis of a, and least for one parallel to the axis of z, or that the axes 
of x, y, z are the axes of greatest, mean, and least elasticity. 

5. If we suppose a range of particles, in a plane parallel to that 
of yz, to be moved upwards, or parallel to the axis of z, the 
resultant of the attraction of all the others will tend to pull it back 
along the same direction with a force proportional to the displace- 
ment: the same thing will happen if it receive a displacement pa- 
rallel to the axis of y, but for the same displacement the force of 
restitution will be greater than before, as we have supposed the 
elasticity for displacements parallel to y to be greater than for those 
parallel to z. Let a, 6%, c® be the coefficients of elasticity, so that 
in the first case the force tending to bring back the range of parti- 
cles is c? x displacement, and in the other 62 x displacement. If the 
influence of one range of particles does not extend beyond the 
nearest range, the force of each of the adjacent ranges will pull the 
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: ‘ , c* : 
displaced range back with a force = > x displacement. 


Let us now consider a succession of ranges to be displaced by 
different quantities, all parallel to the axis of z, the ranges being all 
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at equal intervals. Let the range at a distance a from the origin 
be displaced by a quantity z,, that at a distance x + / by a quan- 
tity z,,,,and soon. The mage at w will be pulled up by the 
range ata +h witha force 4c? (z,,, Z)j= 4 c? Az,, and ‘down 
by the range at «—h witha “force. n Az, »: thus from these two 
it will be pulled up with a force equal to their difference, or 
2 ( tae 2A2> 
$c? (Az, - Az ) = 4 A<z2 
and if M be the mass of the range at a, we have for its motion 
d?z, ce A? 
dt* 2M “anh* 

6. The general solution of this equation of mixed partial differ- 
ences has been given in Article II. of our last Number; and there it 
is shown, that the velocity of propagation depends on the length of 
the wave, and thus the phenomenon of dispersion is accounted for. 
In the present case we shall take the approximate solution, where / 
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a*2 
1c Worv eme 2a ; * 2 for 2+ _¢e oti ee Ig 
is very small. By putting a h? for A®z _,, our equation becomes 


rnmh rah? 


d az = h? dz, 
dt 2M da?- 

This is the well-known equation, which shows that a wave will 
be propagated along the axis of 2 with a velocity Tae . Ifthe 
vibrations take place parallel to the axis of y, there will be a wave 
propagated with a velocity 6 ——. Fresnel states, that he has 
satisfied himself, by experiment, that the velocity of a plane wave 
depends solely on the direction of the constituent vibrations, and 
not on the direction of the wave. It would appear from this, that 
the factors of 6 and ¢ in the above expressions must be constant for 
different positions of the wave; and it results, that plane waves, 
consisting of vibrations parallel to the axes of 2, y, or z, will be 
propagated perpendicularly to themselves with velocities which are 
proportional to a, b,c, or to the square roots of the coefficients 
of the elastic force for the disturbance of a single range. 

7. If the range of particles in the plane of yz receive a dis- 
turbance in that plane, but not parallel to one of the axes, we must 
resolve this displacement into two, one parallel to each axis, and 
there will result fvo waves, propagated with velocities b and e re- 
spectively along the axis of 2, and corresponding waves propagated 
in a negative direction. 

And this result follows, whether the displacement takes place 
along a line, as in polarized light, or in some more complicated 
manner, as in common light; and thus any kind of disturbance in 
one of the coordinate planes will produce two waves, polarized at 
right angles to each other, and of different refrangibility. 
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8. We have now to consider the general case of a plane wave 
inclined at any angles to the coordinate planes. A vibration in this 
plane will produce a force of restitution, which will have a different 
direction from the displacement. Resolve this force into two, one 
along the direction of displacement, the other perpendicular to it ; 
if this last portion be also perpendicular to the front of the wave, it 
may be neglected, and a wave will be propagated with a velocity 
proportional to the square root of the elasticity resolved along the 
direction of vibration ; if it be not perpendicular to the front of the 
wave, it cannot be neglected. We have, therefore, first to inquire 
whether this is the case for any direction of vibration in the plane ; 
and this leads us into the subject which was discussed in our last 
Number. 

A. S. 


Il.—_CIRCULAR SECTIONS IN SURFACES OF THE 
SECOND ORDER. 


In determining the circular sections in the Surface of Elasticity, 
Fresnel has made use of a method which is very readily applicable 
to surfaces of the second order; and as it has not yet been intro- 
duced into any work on Analytical Geometry, it may be useful to 
insert it here. 
Taking first the surfaces which have a centre, let their equation 

be 

et Sh 6 ee ee eee 
and let this be cut by a plane 

2 AON canta sh corber sade aseun (2), 


which we suppose to pass dea the centre, as <a sections made 
by parallel planes are similar. Let this plane also cut the sphere 
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a PP 8 nn cavsnesncsscvnenescisas Cb 
Now, as m, ”, r are indeterminate, we can so assume the position 

of the plane and the magnitude of the sphere, that the circular see- 
tion of the surface (1), if it exist, shall coincide with the section of 
the sphere ; and if these coincide, the equations to their projections 
or the plane of xy must be identical, which gives us conditions for 
determining mand x. Substituting for z in (1) and (3) its value 
from (2), we get 

(P 4+ P’m?) 2? + (P’ + P’n?) y? + 2P"mnay = H...... (4), 

(1 + m*)2? 4+ (1 4 m2) y? + Qmnzy =7"..,.... (5). 
Comparing be term separately, those involving xy will coincide 
if either 
H 
=p 
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Taking the first condition and comparing the other terms, 


H — | r+rr - 1 + n° 
P =—PT", ant 7 =— 2 - 
which gives P’ + P’n? = P(1 +4 2°), 


P— Pp 
and 2 = + Van = 


If we suppose = 0, we find in the same manner 


ms 
a= + 
= | P a P’ ° 


The third condition leads to no result, and therefore is not to be 
considered. 

In the ellipsoid, P, P’, P” are all positive, and 

Parc i. 

This shows that the value of 2 is impossible, and that of m possible ; 
therefore there are two directions arising from the doubtful sign 
in which the ellipsoid may be cut in circular sections, determined 
by the equation to the cutting plane, 


pe 


P - 
s=F Pp’ — Pp av 


In the hyperboloid of one sheet P” is negative, and the value 
of n is possible and m impossible. In the hyperboloid of two 
sheets P’ and P” are both negative, and P” < P’, m is possible 
and 2 impossible. It is true, that for a plane passing through the 
centre the section is impossible, but a plane drawn parallel to this 
at a sufficient distance from the centre will cut the surface in a 
circle. 

The equation to the surfaces without a centre is 

py? + p2? = pp'e. 
Let this be cut by a plane 
x= mz + ny, 
which also cuts the sphere 
zr 4 y? 4 22 — Qe. 

The equations to the projections of the sections on the plane of 

zy are 
py a2 gaan fos <= (Qj ' 
y + pz mppz — nppy = Y, 

(1 + n?) y? + (1 + m?) 2? 4+ Qmnzy — 2mrz — Znry = 0. 

In order that these may coincide, the term involving zy must 
vanish, which will be the case if m=O or n=0. 


If m = 0, then 1 + n° =? and n = + \/® —. 
P Pp 


i ae 
If n = 0, then 1 +m = P and m= + ea ; 
Pp P 
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In the elliptic paraboloid p and p’ are both positive, and accord- 
ing as p’ is greater or less than p, the first or second is to be taken, 
the other becoming impossible. In either case there are two series 
of circular sections corresponding to the positive and negative sign. 

In the hyperbolic paraboloid p or p’ is negative, so that there are 
no sections in which it is cut ina circle. This would appear also 
from the nature of the surface, as it can never be cut by a plane in 
a closed curve. 

The same method may be applied to the oblique cone, so as to 
determine the sub-contrary sections. 

By determining the cireular sections in this manner, it is seen at 
once that any two belonging to different series are situate on the 
same sphere. 


F. 


III.—PROPOSITIONS IN THE THEORY OF NUMBERS. 


THE proposition, “The continued product of any m consecutive 
integers is divisible by 1.2.3...m,” has been strictly proved by 
induction ; but this mode of proof only shows, if we may so speak, 
that the proposition must be true, without showing why it is true. 

We have seen it reasoned, that since the product of any m con- 
secutive integers is divisible by each of the natural numbers up to 
m separately, it is divisible by their product: but a little consider- 
ation will show that this does not follow. 

The problem from which the above proposition is here deduced, 
has been proposed and solved before, but not, that we are aware, 
in the same form. 

Prop. 1. If p be a prime number, the index of the highest 
power of p that is a divisor of the continued product 1.2.3...M, is 
M — S (M) 


p-! 


expressed in the scale, whose radix is p. 


where S(M) represents the sum of the digits of M 


We shall adopt the notation used by several Continental mathe- 
maticians, M! for 1.2.3...M. 


Let the quotient of M by p be M,, and the remainder a 
M = Mp } a> 
and the multiples of p in the series 1, 2, 3, &e. M, are 


0? then 


Pp, 2p, 3p, &e. Mip; 
wherefore, if K denote the product of the remaining numbers in 
the former series, 


M!=K x p2 ».3p...Mp, 
= Kp™1.M! 
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Let 7 be the index of the greatest power of p in M!}, 
un 
then from the last equation 
r= M, +f, 
For the same reason, if 


M,=M,p+a,, and M,= M,p + a,, &e. 


and 7, 7, &e. be corresponding quantities for M,!, M,!, &e. 
71, =M,+ 1%, 
f= M, } %,, 


» 
2 3 


- M 


m= m t m 


r =0, 


m 


7 


M,, being that quotient less than p at which we must finally arrive. 
Adding all these equations, 
r = M, + M, + M, + ...... + M 
But, adding the equations 
M = Mp + 
M, M,p + 4,; 
M, = M,p + a,, 
M,_, = M,p + @ 
M. =a 
and observing that a,, a,;, &c. a@,, are the successive digits of M ex- 
pressed in the scale of p, we have 
M + M, + M,+...+ M,,= (M, + M,+...4 M,,) p+ S(M); 
_M—S(M) 


whence M, + M,+ ..+M,, = eae] ? 


m=~\? 


m? 


and the proposition is proved. 
The following proposition is also useful. 

Prop. 2. If M and N be two numbers, of which M is the 
greater, and S(M), S(N) the sums of their digits to any radix 7, 
S (M —N) is not less than S(M) — S(N). 

For let M =a,7™ 4+ 4,,_, 7" + 0. + Gar” + a, 7?" + «0 

and N = br tb yrt + «ns 
then 
M—N=a,7" +a 


pV (a, —b,) 9" + (dn_y— bn "| + 


m~\ 


Now, as long as the digits of N are not greater than the corres- 
ponding digits of M, the digits of M — N will be the excesses of 
those of M over those of N. But suppose 6, greater than a,, then 
we form the digits of M—N by writing that part of the difference 
thus, 


(4,4, — 5,4, — 1) rt + (7 + a, — b,) 17, 
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and the sum of these two digits is 





Sen, b.., ta,—b+r—1; 
hence S(M — N) may be greater, but cannot be less, than 


S(M) — S(N). 


by 1.2.3.. m. 
N-—S(N) 


p-l 


By Prop. 1, the index of p in N! is 
N+m—S(N+ ) 
p-l 

(N + 1) (N + 2)...0N + ™) on 
is the difference of these quantities, or 


m—S(N 4+ m) + S(N) 


, and in (N-+-m)!, 
; hence the index of p in 


_(N + m)! 
N! f 





p— i 
m — S(m) 
p 


The index of p in m! is ; and, by Prop. 2, making 


therefore m — S(m) is not > m—S(N +m) + S(N); 
and the index of p in m! is not greater than the index of p 
(N + 1)(N + 2)...(N + m). Hence, if 

(N 4+ 1) (N + 2)...(N 4 m) = 2%.9°.5°..., 


OR KY 


and 1.2.3 ...m == 22.38.57 ..., 


—e 


n 


~ 


ais not > a, B not > b, ynot >, &e. Consequently 


(N 4+ 1) (N 4+ 2)... (N +4 m) is divisible by 1.2.3 ... m. 


S. S. G. 


IV._NOTES ON FOURIER’S HEAT. 


Tue method employed by Fourier to integrate the partial differen- 
tial equations which occur in the Theory of Heat, is to assume 
some simple form of a singular solution, and afterwards to extend it 
so as to include all the circumstances of the problem. It is in ef- 
fecting this that he has displayed the great resources of his analysis, 
and imparted so great an interest to his work by the variety and 
ingenuity of his methods. Indeed there is a freshness and origi- 
nality in the writings of Fourier which make them in no ordinary 
degree arrest the attention of the reader. But however much we 
may admire the means by which Fourier has overcome the diffi- 
culties of the problems he had to deal with, yet it seems more 


Prop. 3. The product of any m consecutive integers is divisible 
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agreeable to the usual style of mathematical investigation to deduce 
a result by limiting the general solution by means of the conditions 
of the problem, than by extending a particular case. 

That this may be sometimes done with even more readiness than 
by Fourier’s method, will be seen by the following solution of a 
problem given in p. 161 of the Théorie de la Chaleur. We may 
remark, that there is in general no difficulty in the solution of the 
partial differential equations, but only in the proper determination 
of the arbitrary functions in the solution, so as to suit the con- 
ditions of the problem. 

If a rectangular plate, bounded by two infinite parallel edges, 
have one of its extremities kept at a constant temperature 1, while 
the infinite edges perpendicular to the heated edge are retained at 
a constant temperature 0, the equation from which the temperature 
is to be determined is 

dv d*v 0 

a + er ese (1), 
where v is the temperature at the point X; Ys the origin being at the 
middle point of the heated edge, the axis of x bise cting the plate, 
and the axis of y parallel to the heated edge. For the sake of 
shortness Fourier represents the breadth of the plate by =. 

The solution of the equation (1) by the method of the separation 
of the symbols of operation from those of quantity, is 


v = Cos (y 5) @ (x) + sin (y 4 ) ¥ LS eee (2), 
o (x) and (a) being arbitrary functions of x And it may also 
be put under the form 
vo=FapyV¥—1)4f(x—yV—]}), 
where F (x) = } $9 (x) + Y(a)} and f(x) = } $¢(a) — Y(a}. 
Now, on looking at the circumstances of the problem, it will be 
seen that it must be subject to the following conditions : 


1st, » must be symmetrical with regard to y and — y. 
T v 

L oemé@ehmam — , Whatever @ may he 

2, y =0 when y= 3 5 hatever x may be. 


3, » = 1 when x = 0, whatever y may be. 
4‘, » must be very small when 2 is very large. 


From the first condition we must have (a2) = 0, as otherwise 
the second term would change its sign when —y is put for y. 
Hence we have only 


d 
v = COS (y 5) oN 2 eee (3). 


By the second condition, putting - > for y in equation (3), we have 


O = cos (5 2) 54 errs 


















106 





Notes on Fourier’s Heat. 


Now this is in fact a linear differential equation with constant 
coefficients, and of an infinite order. By the principles laid down 
in Art. V. of our first Number, we can integrate this equation if we 


. , 7 
know the roots of the equation cos (5 :) = 0. Now these are 


= 1, = 3, ca 5, &e. 
being in number infinite. Hence the solution of (4) is 
\ 2 3 be -32 | re. 
PC) nk i i i i te | 
W Hf woe W or 70 
+ C\é€ + C\« 74 Cie? + &e. f 
the number of terms and arbitrary constants being infinite. By 
the fourth condition it appears that the second line of (5) must 
disappear, as otherwise v would be very large when 2 is very large. 
Hence we must have 
Cc’ = @, e. = QO, C. =0Q0, &e. 
and equation (5) is reduced to 


o (x) = a + Cae" + Co? Ae 80.0.0 (6), 
and v becomes 
1 : ro 
v = COS (y <) (Cye* + Cye* 4 Cie? 4 &e.)...... (7). 
ax. ad ‘ 


By the third condition v =1 when «= 0. If then we expand the 
symbol of operation in (7), operate on each term separately, and 
then make 2 = 0, we shall find 


1 = C, cos y + C, cos 3y + C, cos Sy + &e......(8), 


“ie 
and + 35; 


Tv 


where y is contained between the limits — 5 
In order to determine the arbitrary constants, we shall follow 
Fourier’s method of definite integrals. If we multiply both sides 


sa) ‘ — T Us 
of (8) by cosy dy, and integrate between the limits + 9 and — 5° 
all the terms except the first will disappear, as they can each be 
decomposed into the cosines of even multiples of y, which, on 
integration, vanish at both limits. Hence we have 


+ 33 + hn + 30 
fx cos y = cxf cos? y = Sav (1 + cos 2y), 
- 47 —4hr — hr 
whence we find C, = : . 
Tv 
In a similar manner we should find 
1 4 1 4 1 4 
C.2---, C=o-— : a 
° 3m F 5 C, 1 


and so on. 
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Substituting these values in equation (7), it becomes 
T d ee ie, i 
v = cos ly —)} (eo * — we tae Pt — gt &e. ). 
4 (y “.) ( 3 5 7 + ) 
Now if we expand the sign of operation, and apply it to such a 
term as e~"*, we shall find that it becomes cos nye~"*. Hence the 
expression for v becomes 
wv ‘ek od 1 War ¢ 1 52 4 " 
i vse "cosy — 5€™ cos 8y +75 * cos Sy — &e. 
which is one form of the solution which Fourier gives. It may 
easily be reduced to a more simple form. For if we substitute for 


the cosines their exponential values, we have 
ern sat metas BET = 
fe (ay We ans 5 £ 3(arvy 1) + t ge 9(@7Y . = 


b] }~—-4 _— an _— 
L4 e (aty 1) tee 4 g-Xrty=1) -}- cs Pam atyN—1)_ &e. 


3] 


zt 
which, by Gregorie’s series, become 

T = - 
a v = tan ly (x yV~) tan le (aty™ —1) 


<g-yV =} - = 

pe yD + (atyV~1) 

== tan -— ay : 
Le > 


2 cosy 
nk aa 2 
= tan & t.), 


E 
which is the simplest form that the expression can assume. 


D.F. G. 


V.—DEMONSTRATIONS OF SOME PROPERTIES OF 
SURFACES OF THE SECOND ORDER. 


To the Editor of the Cambridge Mathematical Journal. 


S1r,—I observe the author of a Paper in your first Number 
says, he is not aware that any person has made use of the 
symmetrical form of the equations to a straight line: he will find 
them employed in Cauchy’s “ Legons sur les Applications du Calcul 
Infinitésimal.” 

Perhaps the following solutions of two common problems may be 
acceptable to some of your readers. 

1. To find the locus of the bisections of parallel chords in the 
surface whose equation is 

Aa? 4 By? + C2? =k. 
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Suppose any one of the chords meets the surface in the points 
whose coordinates are By Yo Zp Ly Yyp %3 then if a, (3, y are the 
angles which the chord makes with the three axes, we have 


=b. 


cos a cos/f2 ss COS 


zt, — 2, = a 





Let x, y, z be the coordinates of the middle point of the chord ; then 


x a. = 9 \ Yo + Y, zy + zy 


—_———_—>", — > 


2 ‘ ) - 2 


or X + YW + "7 __ afi) + ad 


x y z 


xr — 





; 
therefore, combining these equations with (1), we have 
a il tne ~tawe 
vy — 2% Yo H _ “1 


ec SE ( 


reosa ~ ycosB ~~ zecosy 


Now, from the equation to the surface we have 
eo 2 2 oo ae ar 
A (x, x) + B (Yo — Wo) + C (z,? — 2,2) = 0, 
which, combined with (2), gives immediately 
Ax cos a + By cos B 4 Cz cos y = 0, 
which is the equation to the locus required. 

2. To find the locus of the intersection of tangent planes drawn 
parallel to any system of conjugate diameters in the same surface. 

I shall first observe, that the projection of the diagonal of a 
parallelogram on any line whatever, is evidently equal to the sum 
of the projections of any two contiguous sides on the same line. 
Whence it easily follows, that the projection of the diagonal of a 
parallelopiped is equal to the sum of the projections of any three of 
its edges, of which no two are parallel to each other. 

Hence, if  n, ¢ are three straight lines drawn from the origin to 

Ways) j ; Pp > > y ee 6 rdlinates are xv ~ 
three points P,, I » P,, whose coordinate S are By Yu 2% Ty Y, 2s 
La Yo» Z, and if through the points P,, P,, P, we draw planes 
parallel respectively to the planes of Z, &Z, &, these three planes 
will intersect in some point P, and a line drawn from the origin to 
P will be the diagonal of a parallelopiped, of which ¢, n, ¢ are three 
edges. And if we call x, y, z the coordinates of P, and project the 
diagonal and the three lines & », ¢ on the three axes successively, we 
have evidently, by the principle above mentioned, 

r=%+27,4+2, 
Y=Ytnty, [ove (3). 
P= atte 
Now suppose & », ¢ are three semiconjugate diameters of the 
surface proposed, we have of course 
. 2 eo . 
Az,x, + Byy, + Cz,z, = 0, 
Arr, + Byy, +¢ 27, = %, 
Ax«v,+ By, + Cz, 0; 


” 
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and if we multiply these equations by 2, and add them to the 
three following, viz. 
~ 2 2 ~~ 
Ax, + By, + Cz, =h, 
2 2 = 
Az,? + By? + Cz? = &, 
Az? + By,? te Cz, — hk, 
we obtain immediately 
= Vee ee ae ey | 
A (at) +2,+2,)? + By ty, +Y,)? + C (2, +2, +2,) = 3h, 
or by (3) 
Aa? + By? + Cz? = 8h, 
which is evidently the equation to the locus required, 

I do not know whether this solution is new, but I do not re- 
collect to have seen it before. 

The equations (3) lead immediately to the common formule for 
the transformation of conndinaton; for if we call ay By yo ay By Yp 
tg) Bay Y» the angles which é, 7, ¢ make with the three axes respec- 
tively, we have x= écosay y, = §cos By z= cos yo, &c. and 
the equations (3) become 

x = Ecosa, +7 Cos a, + ¢ cosa ty 
y = & cos 3, + 7 cos 3, + ¢ cos | 
z= cosy, + n cos y, + f cos ‘~ 
where & n, { are evidently the coordinates of P, referred to lines 
drawn from the origin through Py, P,, P, as axes. 
I am, Sir, &c. M. N. N. 
Oxford, Jan. 31, 1838. 


VI.—ON GENERAL DIFFERENTIATION.—No. II. 


In our first Number we investigated the general differential coeffi- 

cients of e"* and 2". We shall now proceed to find those of other 

simple functions, and to give some examples of the application of 
the theory. 

d+z” 

1. We found that our formule gave an infinite value for aa 





a? 
when x is not, and x—a is, a positive integer. But by using the 
complementary function, a finite value of another form may be 


obtained, in the same manner as the value a [FT may be derived 


from that of fx"dz. 


Let a — in formula (N), and x — a in formula (O), (p. 19), be 
assumed equal to r+/3, when 7 is an integer, then the second 
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utp 


member of each will be a multiple of » Which is 


sin (r + B)= 
infinite when 8 =0, but by taking in a part of the complementary 
a” (a8 — a8) 


sin (r+) x’ 


x , 
#” log - when3=0. Hence in these cases (N) and 
a 


function becomes 





which is a vanishing fraction whose 


Sad’ al 





value is ~ 


(O) become respectively 


ee l z 
i= =(-1)™ ——— #2" log — ...(Q), 
dx~" x" ( ) rd +a—n)(n) a (Q) 
d“x" sin r I'(1 + 2) " x 
= (— 1)" —— . ———_~., 2 log - ... (R). 
dax ( ) us r(l +2—a) “ 8 a 7 
“loo x 
2. The value of : i may be readily derived from the 
da* 
formule for the differential coefficients of powers of x. For since 
d — ] 
or z = - 
dx ° a’ 
d* ee a*-' 1 
dee 8“ = Gai x" 
Hence, if a be positive, 
d* — a—l I r ) Y 
dx" log « = (— 1)*-'T(a) > eka eieanied (S), 
° iia 7 : - 
, ee Gee | ——_— 2, by (N)...... r). 
ani dx~" log z 9 snaml'(1 +a) — 7) ( ) 
3. If it be known that 
d+ 
dx ¢ (r) tot y (*), 


it is evident that 


— (mx + a) = m*) (mx + a); 


hence the general differential coefficients of all rational functions of 
x can be found, by preparing them as for integration. Also, if the 
differential coefficient of a function to any integral index be a 
rational function, its general differential coefficient may be found. 
As an example, let 





y = tan™' 2, 
dy ] ( l 1 
ea / —" + a ); 
sane dx ] + 2 7\1 Ly — 12 1—V—I12 
therefore, using the formula (P), 
P(— Df (7—1)° 
(14% — 12)" 


a’ phe Ff 4}. 
(1 ~Y—12)"J 


-~ di 
dx" = P( a) 





eee 











ry 


se 


id 
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Substituting for w its value tan y, observing that 
¥—1 = cos (2r + } 4) 3 +V—1 sin (2r + $)7, 


and reducing a ‘ Moivre’s theorem, we obtain 


a= - Fr + _— y) cos a sy — (2r + h) a n. 
P(- i d= ie 

The ws >» of — s the - la " 
The value of P(—a) is that of a. oc which i 


(—1)*? P(a) 


when a is positive, and 


. T 
(—1) sin (— ar ).I'(1-— a) 


when a is negative. The formula fails when a is a negative integer. 





4. Let us ascertain in what cases the differential coefficient of a 
constant quantity is not zero. 
Since C = Ca, 
d* C P (0 e C i 
oi OO ig a cat 
daz P(—a) P(—a) 
which will not be zero if P(— a) be not infinite, which is only 


when —ais a positiv e intege r, oraa negative integer, that is, in 
the case of common integration. 


5. The differential coefficient of infinity may be finite. For 








2 = CP(m) x", if m be not a positive integer: therefore 
q"= , P(m) P(n) pone 
dx* ~~ P( (n - —a) ” 


which will be finite if x be a positive integer, and x — a not. 


6. We proceed to find the general differential coefficients of 
cos # and sin 2. 


Since cos + “W— 1 sin a = ev! * 


dat 
= feos (2r + 3)a + Y —1 sin vl + 3) 724(cos a + 7 — 1 sin x) 
= cos $(2r + 4) am + at + ¥—T1 sin §(2r + 4) az + 2. 


In like manner, 


da? cos x PY i@ * sin vc =(V—] se ips 


d* cos x PY om 7 d* sin x 


dat - da* 
= cos $(2r' + L)ax + 23 —V—} sin $(2r' + 4) ax + 23, 
r being some integer, not necessarily the same as 7. Adding and 
subtracting, 


~ 
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d* cos x 
da 
= feos (r — 1’) aw 4 v7 —1sin (r — r) ane 
cos f(r + 1! + By aw 4 2b, U 
d* sin x p sevens) 
dx 
‘ / s ’ 
= feos (r — r) am + V — I sin (r — 7’) an} 
sin {7 +) + 4) ax 4 23,7 
Since r and 7’ are both arbitrary, there is no relation between 
r +7’ and r —7’, except that they must be both odd or both even. 


If a be a rational fraction whose denominator, when in its lowest 
terms, is 2, 7+ 7’ must go through all the 2n values, 


0, 1, 2 3, 00... Zn — ], 





before the variety of values of the trigonometrical expression into 
which it enters can be exhausted; and the same is true of r — 7’. 
Now any of the » odd values of 7 4 7° may be combined with 
any of the 2 odd values of 7 — 7’, and thus n® different values 
of the total expression are found. The even values of r 4+ 7" and 
x —r produce as many more different values, hence the fractional 
differential coefficient of each of the quantities cos 2 and sin a has 
2n* different values.* 

For example, the differential coefficients to the index 4 of cosa 
and sin a have 2 x 2° or 8, values. Those of the former will be 


found to be 
7 3r 
+ 00s (7 +2), £7=1 00s (= +2), 


53 wks 
+ cos > + r), + /— 1 cos (= + r). 


The reader may observe that if differentiation to the index 4 be 
performed twice in succession upon cos a by any one of these 
formule, provided the same be used both times, the result will 
always be —sin 2, as it should be. 


7. The same multiplicity of values attends the differential co- 
efficients of other trigonometrical functions. For instance, let us 
consider the functions «”* cos ma and ¢”* sin nx. It will be con- 
venient to asume m = a cos 0, and xm =asin@. We may then 
change az into z, and they become 


e* 59 eos (x sin 0), ¢* ©°S9 sin (x sin 6), 


which we shall call © and 0’ respectively, as Mr. R. Murphy has 


* M. Liouville did not employ the direct process which we have, and 
consequently did not discover that the differential coefficients had more than 


2n values. 
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done in the Camb. Phil. Trans. vol. v, where he has proved some 
properties of them. We have then 
6 + V¥—1 0! = &* 9 feos (a sin 0) + VW —1 sin (w sin 092 
= gt (cos 8+¥=1sin 9); therefore 
d 
dixt 
= fcos a (2rr + 0) + V¥—1 sina (Qrr + 0)3 
e* © 8 Seos (a sin 8) + 7 —1 sin (a sin 0)3 
= eT 89 Scos (a sin 0 + Zarz + af) + sin (asin 0 4 2arz + a0). 


In like manner 


(0 4 V/—1 0’) — (cos 0 4 eg sin 0)" g& (cos 6+ /—TI sin 6 


[a anne 
ie (0 —-V—16') = 


e7 89 Sos (x sin 6 -+ Qar'r + af) — sin (asin 4 Zar'm 4 al). 
By addition and subtraction 
a0 
dx 
Scos (r — r) ar + V— 1 sin (r — 7’) an} 
e7 © 8 eos fa sin 8 + (r + r) ax + ab}, 
d*0’ ae \ 
dat | 
Scos (r — r') am + V— Lsin(r — 7’) an} | 
e7 © Osin Se sin 8 + (7 + 1’) am + al}, , 
As before, the number of values of each of these expressions is 


twice the square of the denominator of a, supposed a rational frae- 
tion in its lowest terms. 





The necessity of supposing r and 7’ different, may be proved as 
follows. If 0 be changed into — 0, 0 remains the same and 0’ is 
altered only in sign: no new value of the differential coefficient 
ought therefore to be found by making this change: but that would 
be the case if 7 and 7’ had been supposed equal: therefore that 
supposition would not give all the values of the differential coeffi- 
cients. A similar remark applies to the last article. 


8. We now come to another division of our subject, in which we 
shall prove certain formule discovered by M. Liouville, and ap- 
plied by him to the solution of a variety of problems. The first 
of these is 


@ d-" 
fe (x + a) a"™"' da = (— 1)" T(n) pm $ (2) .0.200006(W) 
0 


This may readily be proved by the theory of generating functions, 


Q 
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Let G@(wx) denote the generating function of g(x), or that 
function of ¢ in the expansion of which the coefficient of ¢ is 
(x); then we have 


G o(x + a) = €°G g(x), 
G/o(« + a) ada = ft a"! daG 9(x). 
0 0 






























0 : 
Let @ = 2, then a= porns and the second member of the 

or 

fe} 


preceding equation becomes 


['(n) (log t)-" G p(x); 


] dr 
but (log ) ‘'Go(«#) =G ae * (x): 


o "il 
therefore G/@o(# + a) a""'da = (— 1)" T'(n) Ga iam t 4 (2), 


2 d- 
and J ) (x $+ a) a” ‘da = (— ] , by (n) dae rH) (2x). 
0 


It may be remarked, that if the form of ¢(2) be such that the 
first side is infinite, the second may be made to agree with it by the 
aid of the complimentary function, if it do not without ; so that the 
formula is true for all forms of p(x). 


9. An obvious use of this formula is to find the form of ¢ (x) 
when the value of the definite integral is known, as in the following 
problem. 

Each particle M of a line AB, infinite in both directions, exerts 
upon a particle P without it a force perpendicular to the plane 
ABP, and proportional to sin PMA, and to an unknown function 
F(r) of the distance PM or 7. The total action of the line upon P 
is a known function f(y) of the perpendicular distance y of P from 
AB. It is required to find the function F (7). 

This is an important physical problem, for if AB represent a 
voltaic current, and P a pole of a magnet, the action of M on P 
is such as has been described; and the total action is found by 
experiment to be inversely as y. 

Let C be the foot of the perpendicular from P on AB, and let 
CM =z, then r* = y* + z*, and sin PMA = y whence the 

rT 


whole action of AB is 
+e 
ry FV (y? + ar. PVG 4D) et?) a 
ye 2 ° 2y VP 42 ie 


“== 0 


QV 
| 


which is equal to f(y) by hypothesis. 
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), 
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F (7) 


Assume — ss »(7*), therefore 
er ae Sy) 
2/o(y +. 2°) dz = ; 
0 y 
Let y? =a, z* =a, then the first member becomes 


/ ‘@ (@ 4 a) am da, 
which by (W) is equal “ 


’ A= 
ee r(3) (;.) g(x); 
fi 


— 1\-3 
therefore VW — 1 r (4) (=) ¢(x) = - 


. ia l SV 2) 
and ¢(«) = V_1Tr (4) (3) yu 


Supposing that the differentiation is performed, and 7* is substi- 
tuted for x, F(r), which is — to 7 (7°), is known. 


If f(y) = sa fv =$ == and by formula (E), p. 14, we find 


ois = tax™ y 
; a 
therefore F(r) = 5 


It is easy to verify this result by ordinary integration, and thus 
to confirm the truth of the principles of the new calculus. 


10. It will be convenient to possess a formula in which the 
limits of the integral shall be 0 and 1. Such a one may de ob- 
tained from — by changing the independent variable. 


Let a+a=~—; then when a=0, 0=1, and’ when a=o, 0=0: 


1 — ad0 . : p 
also, a = ns % d= _— , and the formula is changed into 
1 
x a" n » 
[o (5) gen 0-80 = IP) TS 6) 


0 
Let at g(a) = (- ): therefore 


fC :) (1—0)""! do=(—1)" F(n) « alan J (;) p09. 


0 
The following modification of the formula is often convenient. 


] 
Assume 0 = fa and # = =, then we have 
Y 


/ WB) (y—)"— dB =(—1)* Pm) yd Sy" WY (y 2 (a. -) or (9° 
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The differentiation indicated in the second member must be per- 





, : ] , : 
formed with respect to — as the independent variable. 
y 







11. Asan example of the use of the last formula, we select the 
following problem, which includes, as a particular case, that of find- 
ing the tantochronous curve. 





To find the curve such that the time occupied by the descent of 
a particle sliding along it by the action of gravity to a given point 
from a given vertical height above it, shall be a given function of 
that height. 




















Let x be measured vertically. Then the known expression for 


the time of descent from a height / to the origin, is 
h 


1 ds i 
a ——— (2 * 
les (h nr x) dx 


0 
but by hypothesis this is equal to the known function f(/). 
U: , : , . 
Putting i = y(x), and comparing the integral with that in 
dx 
formula (Y), we find 


VITA a? H W(AR (a.3)° = 9g f (h), 


whence, observing that 1 (4) = 7, 
a oF ee 1\~4 
VA) =— V1 4/8 a yap (a.;) z, 
T fis h 
By changing h into 2, (2) is known, and thus we have a diffe- 
rential equation to the curve. 
Suppose, as a particular case, that S(h) = ch"; then, after 


changing 4 into - for convenience of operation, we have 
2 


l — 2 a2/fd\i “5 
¥(:): ike Vee i 


id / ¥ lr (n + 1) wr 
— € TD (n+4)° 


by (F), if 2 be positive. Therefore the equation to the curve is 


ds —_— 29 T (x +1 ) yn-3 
de ~° Vox P(ath)* ” 


If » = 0, or the time be independent of the height, we have 


ds 29 _, 
—_—- =C — xX 8, 
dx x 


which belongs to a eycloid 


i 











per- 


the 
nd- 
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of 
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If » = 4, or the time vary as the square root of the height, 


ds Dae 9 
dx °'V/ 9° 
which belongs to a straight line. 


If x = 1, or the time vary as the height, 
d 27% 
is _ 2V% ai: 


=c 


dx T 


2V%_ it 
T 


putting ec = Ta , and solving the equation, we obtain 


2(«#— aye 
ee es 
at 
We might give many more instances of the use of these formule, 
but our limits will not allow us: and therefore we recommend to 
such of our readers as are particularly interested in this subject, 


the original Memoirs above referred to, which are more complete 
in the applications of the Calculus than in its principles. 


UE 





VII.—ON THE TANGENT PLANE OF THE ELLIPSOID.* 


In an article, which appeared in the first Number of this Journal, 
some examples were given of the advantage of employing, in cer- 
tain cases, a form of the equation to the tangent of the ellipse, 
which does not involve the coordinates of the point of contact. 
There is an analogous form of the equation to the tangent plane to 
the ellipsoid, viz. 
la + my + nz = VEPa® + m*b* 4. née’, 

1, m, n being the cosines of the angles which a perpendicular on 
the tangent plane makes with the axes. We shall not give a proof 
of this, but refer the reader to Hymers’ Geometry of Three Di- 
mensions, Art. 29. A few examples of its application may per- 
haps be useful. 





1. To find the locus of the intersection of three tangent planes 
to an ellipsoid, which are mutually at right angles. 


Let (dmn), (Um'n’), (l'm'n') be the cosines of the angles which 


* From a Correspondent. 
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perpendiculars on the planes make with the axes. Their equations 
will be 


lx + my + nz = VPa? 4 mb 4 nee, 
Ce + my + nz = VP 4 m2b? + n'*e?, 
Va + my 4 n!z = VIP 4 mb? + nlc. 
Squaring both sides of each of these equations, adding, and 
observing that, since the planes are at right angles, 
P4l?240%=1, m4m24m2=1, n24n24n"=1, 
lm +-Um' 4+ Um" =0, In4lal +n" =0, mn+mi'n' +m'n' =0, 
we have 
a 4 y a 2a + is 4 ce. 
This proposition was first proved by Monge. 
2. To find the locus of intersection of the perpendicular on the 
tangent plane with that plane. 
The equation to the tangent plane being 
le 4 my + nz = VPa® 4 mb? + ve 


the equations to the perpendicular will be 


2 y_32 
Lom Hn 
Multiplying the successive terms on both sides by these equal 
quantities respectively, we have 
ey? +4 2= Var 4 By +4 22%, 
or (a? +y? 4 2°)? = ata? 4 by? +4 cz?, 
3. Three planes, mutually at right angles, touch three spheres 


whose radii are a, a’, a’, respectively. To find the locus of inter- 
section of the planes. 


In this case the equations will be 
le + my + nz = a, 
lx + my 4+ vz = a, 
Vx + m'y + n'2 =a’. 
Squaring, and adding, we shall have 
a 4 y? 4+ 22? — a® 4 a? 4 "2, 
4. Three planes, mutually at right angles, touch three concentric 


ellipsoids, whose principal sections have their foci coincident. To 
find the locus of intersection of these planes. 


‘pee “ . ° ° 
Let abe, abc, a’b’c’, be the semi-axes of the three ellipsoids ; 
then the equations will be 
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lx + my + nz = VPa? - + mb? + nX?, 
te 4+ my 4 vz = VPA? mb? + n%e?, 
Ux + my 4 nz = Via + mY? + ne". 
Squaring both sides, and adding, we have on the first side 
av + y*? + 2°, as before. Also, 
Pa? + mb? + nc? = a? + m?(b* — a?) 4 n* (ce? — a?), 
Lal? 4. m'2B'2 4 n'%'2 — B2 4. 12 (a2 — B2) 4 n'2 (e2 — 8), 
which, since the foci in the principal sections coincide, 
= 6% + 1% (a® — b*) + n* (ce? — B*). 
Similarly, 
Hq !9 4 mlb"? 4. nl'2e!* — cf? 4. 2 (a® — c2) + m"? (b? — c*). 
Therefore the second side becomes, by addition, 
= a? + bY + 
+ a? (2? + I"? — m* — n*) 
+ b? (m* 4 m”* —I* — n*) 
+c (n® + n* — I? — m"). 
But P+ F844 MPs la Py mi 4 nt; Ot + = mit + er’. 
Similarly, 
m? + m'? — U2 4 n'® and n? 4 v2 = I”? + m"?, 
Hence the equation to the locus is 
epy tema + 6% 4+ c% 


Again, since 
em be — a? —_ b'2, 


or a? + 62 = a? + B, 
and 62 + c’2? = b'2 4 @; 
therefore, by addition, 
a+ b24¢c2 =a? 4+ 62 4 &. 
And, in the same way, 
a? + b2 4 2% = a? + B+ 2, 
Consequently the equation may be put into the form 
at fy? fata} f(a +b fe%) (a? 462402) 4(a" 4b" pe"). 
This result was first published in the 19th volume of the 
Annales de Mathématiques. 





If the equation to the surface be 
y* z 
2+ 


= 


2 
= 
the equation to the tangent plane may be exhibited in the form 
L (le + my + nz) + 3 (am? + adn?) = 0. 
As there is no proof of this in Hymers’ Analytical Geometry, 
we subjoin one. 
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Let the equation to the tangent plane be 


This must be identical with the common equation to the tangent 


le 4+ my 4+ nz = 6. 


plane drawn at a point (a'y'2'), viz. 


‘ al 
r + a’ —_ 9 é f =), 


therefore we must have 


a a 
2 "y Q2'z x 
or We We _ 2, 
ax ax 
l 1 ; é 
Ss), orr=-s 
C a’ 7 a 
os wf Qly' 
r) ax’ ao ‘ 
: ma 
whence ek. as ee 
od a 
xs , na 
Similarly, 2’ = — — 
al 
{2 ' 
y zie 
But 2’ = “ 4+ —; 
a 


_ ma, na 


6 
hence —- T— ae + ae 


s 


1 
eins (ma 4+ n*a’). 


Consequently the equation is 


U (le 4+ my 4 nz) + 3 (ma + n2a’) = 0. 


The reader will find no difficulty in proving, by means of this 
equation, that tangent planes, mutually at right angles, will always 


intersect in a plane whose equation is 


and that the locus of intersection of the tangent plane and the per- 
pendicular upon it from the vertex, is a surface, represented by the 
equation 


r+i(aya)=0; 


x(a? + y*? + 2°) +4 (ay + a2*) =0. 





W 





int 
ob 


Ww! 
to 
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VIII.L—TRANSFORMATION FROM RECTANGULAR TO 
POLAR COORDINATES IN DIFFERENTIAL EXPRES- 
SIONS. 


WHEN we have a differential expression involving two independent 
variables, which we wish to transform into one with two other 
independent variables, the method to be pursued is not at. first sight 
obvious. If, for instance, we have the double integral 


SL V dx dy 
where V is a function of 2 and y, and we wish to transform it so as 
to involve rand 0, the four quantities x, y, r, 0, being connected by 
the equations 
2 = fir, 0), y = F(r, 6), 

so that 

a= - dr + - dé, dy = dy dr + dy dd, 

dr do , dr dd 
we cannot multiply the expressions for dx and dy together as they 
stand ; because in the integral y is supposed to be constant when a 
varies, and x to be constant when y varies. We must introduce 
one of these conditions by supposing de = 0 when y varies, or 
dy = 0 when 2 varies. 
Taking the first of these, we have the equations 


dx da 
0o= dr +4 dd 
dr d) 
d. 7 
dy = a dr ey do. 
F dr do 
Eliminating dd between these, we find 
dx F (= dy dx dy\ , 
at = = — — id. 
do “4 \d0 dr dr di] 
From this it follows, that when dy = 0, dr == 0. Hence 
dx 
dx = — dd. 
dé 


Substituting these values in the double integral, it becomes 
.. [dx dy dx dh 
S/ V ( a cs J dr dé, 
do dr dr do 
where V involves only 7 and 0. 

When 2 and y are the rectangular, 7, 6 the polar coordinates of 
a point, x = rsin 0, y = 7 cos 0, 
and S{V dx dy = f{{V rdr do. 

If there be three independent variables the same method is to be 
pursued ; but as we have to eliminate between three equations it 
becomes very complicated, even though we avail ourselves of the 
method of elimination by cross multiplication given in our first 
R 
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Number. In the particular case of the quantities representing the 
coordinates in space of a point, the assumption of a subsidiary 
quantity much facilitates the calculation. For let 
SS[N dx dy dz 
be the function to be transformed, and 
xr=rcos0, y=rsin@ sin g, z =r sin 6 cos ¢. 
Assume p =r sin 6. Then we have 

y=psing, z =p Cos@......... (1), 

p=rsing, 2= cos 6......... (2). 

We can then transform, as we did before, from a, y, z, to 
x, p, o>, and after that from 2, p, @ to 7, 0, 

The simplicity of the method consists in this, that the second 
operation is exactly similar to the first, so that we do not require 
to repeat the calculation, but merely to write down the result, 
substituting 7 for p, and @ forg. Effecting the first operation, 
dy dz) 
do do dp dy | 


{ dy dz 


dx dy dz = dx dp do 


which by equations (1) becomes 
dx dy dz = p dx dp do. 
Again, changing from 2, p, @ to 7, 6, ¢, we find 
dx dy dz = pr dr dd do = r* dr sin 6 dO do. 
Therefore 
SI[N dx dy dz = JS{[N ?* dr sin 6 dé do. 
The same assumption may be usefully applied to transforming 
the expression 
d?V d2v d2V 
dx? * dy? ‘dz 
from rectangular to polar coordinates; for it will be seen, that 


OCT 


' dV d2\ : ’ 
having found —— and ,» the whole expression may be written 
, ) 


dy 
down without further trouble. 
Thus we find 


dV a ii dV L cos @ dV 
dy ? dp p dy : 
ov... d2Vv cos? @ d*V cos? ¢ dV 
=~; = sin? ¢ = —— - 
dy* dp? p> dg? p dp 

2 si ‘08 @ 12V 

mt 2 sin oo q (, C A ” dV 
p* \ dp dy do 


. . ( . . . 
The expression for 72 '8 got by putting 90 — for @ in the 
dz* ; 


above, so that without writing it down 
av  dV_ dV 1 d?V 1 dV 
dy? dz? dp2 p> dg? p ‘dp 
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Putting r for p and @ for ¢, we have the similar expression 
dV  dV_ dV 1 d?V ' 1 dV 


dp “dx drt r de r dr 


ee ia el ; 
And the expression for 7 ur and @ being similar to that for 
ap 


y 


dv . ) 
—— in r and 49, gives 


dy 
Lav _1aV | cota av, 
p dp r dr rr dé 
Adding these three expressions, 
dV d2V d2V 1 /d2V cos 0 dV 
dx? dy? dz - (i sin 0 do) 
l d2V d?V 2dV 


— r 


r?sin?@ dg? v drt r dr 
y J2V 2(7V 
a 1 d (sin2 9 dV ie ] : d*\ ae l d (rv ) 
r? d cos 0 d cos 0 r? sin* 0 dg? r dr 


which is the well known expression. 


IX.—ON THE SOLUTION OF PARTIAL DIFFERENTIAL 
EQUATIONS. 


Tue integration of Partial Differential Equations is much facili- 
tated by the principle which we have developed in our preceding 
Numbers, and it is the more remarkable, that it has been so little 
applied to these equations, as the first step was taken many years 
ago. Fourier, in his 7raité de Chaleur, published in 1822, has 
shown that the series which are obtained in the solution of several 
partial differential equations may be conveniently expressed by the 
separation of the symbols of operation from those of quantity. But 
though he has used this method very frequently, yet he appears to 
have had some unwillingness to give himself up to it entirely as a 
guide in his investigations, as if he were not familiar with the prin- 
ciples on which it is founded. His idea apparently was, that the 
expression which he obtained as solutions might be conveniently 
expressed by separating the symbols of operation, and not that the 
symbolical expressions are the proper solutions of the equations, 
and the series merely the expansion of them. Other French 
writers seem to have avoided carefully entering at all on the track 
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which Fourier opened: Poisson in particular, in a digression on the 
subject of partial differential equations in the second volume of 
his Mécanique, does not put in the symbolical form the solution of 
a very simple equation, which is so given by Fourier. Mr. Great- 
heed, in a paper published in the number of the Philosophical 
Magazine for September 1837, was, we believe, the first to call the 
attention of mathematicians to the utility of this method in the case 
of partial differential equations, but he had not then reduced it to 
its greatest degree of simplicity ; and his paper is chiefly occupied 
with a particular class of equations of the first order with variable 
coefficients, which are not so interesting as many others with con- 
stant coeflicients. We shall here, therefore, proceed to give several 
examples of the application of the principles which we laid down 
in Art. V. of our first Number. 

And, first, we may observe generally, that linear partial differen- 
tial equations between any number of variables with constant coef- 
ficients, are to be treated exactly like ordinary differential equations 
with regard to one of the variables, the symbols of operation of 
the others being treated as constants. If, for instance, we have 
the equation 

dz dz 
a Ix—- = Cc 
dx dy 


This may be put under the form 
d j *,) — 
(« = + -* 2 =, 


ae d 472 -_ 0 
= (« ae oa b=) (¢ +. )s 


and therefore 


: d , 
where we suppose a to be the variable, and - a constant with 
ad, 


: : : d ey 
regard to it. Now the operation (a — + 6 — by the theorem 
dx dy ? 
given in page 25, is equivalent to 


6 ad b i 


C 
Fg - I — 
a Ps a dy [dx e@ dy | 


which being performed, gives 


— =) A. 1 -~ze iv f ? 
-=(5 tart "S(Y); 


JF (y) being an arbitrary function of y, taking the place of the 
constant in ordinary differential equations. 


From this expression we get 


_ 7” 
z= + flay me), 
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as by Taylor’s theorem 


d 
e dz f(a) = f(x +h). 
We might equally well have supposed y to be the variable, and 5 
a constant with regard to it. 
Again, taking the well-known equation for the motion of waves, 
d*z , d*z 


a 


aha Ga ae = O 
dt? dx? ; 


it may be put under the form 


(F a) 7 
ee aa 


and integrating it like the ordinary differential equation 


at cf —at i 
we find z=e« “o(x)+e “ W(x), 


or z = o(@ + at) + W(x — at). 

The equation r — 2a s+ a*t=0 may be put under the form 
d ) 
re 

£3 dy 


the solution of which is by the Theorem in page 25, 


= 0, 


a 


d 
z= WY f2dar.0 
a 
" fro(y) + VY), 
¢(y) and w (y) being arbitrary functions of y arising from the in- 
tegration. Hence, finally, we have 
z= «p(y + ar) + Hy + a) 


The equation 


ar 
€ 


r — at 4+. Zabp + 2a°bq = 0 
is equivalent to 


d d lid d 
— — — —j)|}z—0. 
(z (« dy 2ab) J (3, 7” ra 


Integrating with regard to the first factor, we have 


d d £ (. = —2ab ) 
(5 +? x) zeae Oy) 


Integrating with regard to the remaining factor 
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d d d 
— e "dy fdy oorlaG -b) 4(y) 4 ot dy Wy) 
d d -1 
~ar~ 2ax\~ —5 d } ; 
ae dy (z ) {oa (% — 5 (Y v(y — ax 
Stig E dy Jf oy) + Ly ) 
d 
az onl . 
=e dy gabe by o,(y) + v(y — ax) 
by changing the arbitrary function 
— ,-2abe ,W(y+ar) o,(Y + ax) + viy — ax) 
= 9-42) » (y + ax) + Y(y — az). 
Let us take also the equation 
r— ct = ay. 
Without operating on each factor separately, we may arrive more 


readily at the result by the same means as those employed in 
page 28 of our first Number. For we have 


_ (% con + : 
Gera - o(y + ax) + H(y — ar) 


d2 
~ dx ® (1—« “7 oa =) (xy) + o(y + ax) + (y — ax) 


d2 d-?\-'! 
= i~ eh (3 z) + o(y + ax) + U(y — a2). 
Therefore 


= 4 oly + ax) + U(y — ax), 


the arbitrary functions being added as in the second example. 
A very simple equation being one which occurs in the theory of 

heat is 

dy d*v 

—_—_ s = a — 

dt da?’ 
which is the expression for the rectilinear propagation of heat. 
The solution of this is easily seen to be, if we integrate with regard 
to ¢, 


d2 
. = Pe dx? f(x), 
f(x) being an arbitrary function. If the sign of operation be 
expanded, we shall obtain a series which is the solution derived by 
Poisson from the method of indeterminate coefficients. Laplace 
has deduced from the series an elegant expression for the solution 
under the form of a definite integral ; but it may be more easily 
deduced from the symbolical solution in the following manner. 


Since fran &@ dw = Ax, 
andalso /@, e"°" dw = Wz, 
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we can put the expression for v under the form 


ad? d 
Vr = y em dw gl dzt & (ws ” dr)" (x) 


2 


d 
= [2 dee ovat dz f(x) 
oe 
= FE dw e~* Kz + 2 Vat) 


by Taylor’s theorem; and this is Laplace’s expression. 


Qy 

. dv d% 

This equation — = a, not being homogeneous in the index of 
at 


the operations, admits of two solutions of very different characters. 
The one, which we have found by integrating with regard to ¢, con- 
tains only one arbitrary function of a The other, which may be 
found by integration with regard to x, must contain two arbitrary 
functions of ¢, as the index of operation is of the second degree. If 
we write the equation in the form 


dv tldv 


dv adt — 
and integrate by the method employed in page 28, we find for 
the integral 

v= (= 4. 


is @ l x 


d* . 
2 en ' 
123 a+ a@ 120345 det * ) o( ) 


1 aw d 1 at d2 
+ (1+; 12 atta 1.2.3.4 dé2 + &e.) 40), 

It seems at first sight anomalous, that the same equation should 
have two solutions so different in character: the following is the 
explanation of the difficulty. Since by Maclaurin’s theorem any 
function of a variable may be expressed by meaas of its differential 
coefficients taken with regard to that variable, for the particular 
value 0 of the variable, we know the function if we can determine 
its successive differential coefticients. Now from the equation 
dv ss dv 
ae” dx? 
determine the values of all the differential coefficients with 
regard to ¢ when ¢=0. So that in the resulting expression 
deduced from Maclaurin’s theorem there is only one quantity left 
undetermined, which is the arbitrary function f(x), introduced in the 
dy 1 dv 
dx? ns a dt 
determine from the value of v when x = 0 the values of the alter- 
nate differential coefficients with regard to 2 There must conse- 
quently yi ‘enaemeen ane Yr undetermined quantity, namely, the 


we can, when we know the value of v when ¢=0, 


integration. But from the equation » we can only 


value of * sd when 2 = ‘0, for knowing these two quantities we can 
vr 
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determine the values of all the successive differential coefficients 
with regard to a. 
The equation for determining the vibratory motion of an elastic 
spring is 
dv div 
d@ — dx‘ 


The solution of which is readily seen to be 


]2 
v = COS (« =) F(x) + sin (« ,) f(x). 
l da?/ 


The equation for determining the vibratory motion of elastic 
plates is 
dy | d4v dy d4y 
ar " 4 j — 0 
dt? da’ dx? dy? — dy* F 


which may be put under the form 


d?v ad? i? 
+ (5 +o =0; 
dt? da dy? 


the integral of which is 


d? “a , 
v = cost (5 t =) F(a, y) + sint (+ t 5). f (2, 


In investigating the motion of heat in a ring, we obtain as an 
equation for determining the temperature at any time and point, 

dv ; dy 

a ae hv; 


the solution of which is 


Tee) 
vet x? ) S (a) 
k a 
Pe 
=f tt dx? S(*), 
an expression closely connected with one previously given. 

In the examples which we have given, the coefficients are con- 
stant, but if one of the variables only enters into the coefficients, 
and we integrate with regard to the other, as the one variable is 
unaffected by the sign of operation with regard to the other, it may 
be considered as a constant in the integration. 


A good example of this kind of equation is that which expresses 
the motion of heat in a solid cylinder of infinite length, namely 


dv (= 1 =) 
= @Z —_—_— 9 
dt . dx* F x dx. 
the solution of which is 


7 = 


a ld 
ae (Gants as) F¢n), 
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or, as it may be expressed, 
ld 
;— “ais ( is) ¢(2), 

But this will only be possible when the integration is of the 
first degree, since, as we showed in our first Number, the sym- 
bols of operation become subject to different laws when the va- 
riable itself and the sign of differentiation are both involved. 
And this leads us to the consideration of equations with variable 
coefficients. As in the case of ordinary differential equations, the 
solution of this class is attended with great difficulties, so as to 
become almost impossible for equations of an order higher than 
the first, it is to be supposed that these difficulties are no way 
diminished in the case of partial differential equations. Of these, 
however, when of the first order, Mr. Greatheed has shown, that 
a large class may be solved like ordinary differential equations. 
It is included under the form 


dz a» AZ 

—-+XY—=Pz2+Q, 

dx dy 

where X is a function of 2 only, Y of y only, and P and Q 
of both x and y. This wai” be reduced to the form 


dz 
—_-_= Pa ¢ : 
dx +4 = Pe +Q 
by acl Pim wivele.imti=2 
y a che re > ware ds . = i. of n a at ae 8 
dV a € lange oO le variable; for 1 y [y ’ 1 dy dy 


so that it is only necessary to consider the latter equation. Let 
it be put under the form 


dz 4 (x d . P) — 
dx dy 


and treated as an ordinary linear differential equation between 
z and x. If we integrate by the method of integrating factors, 
J ax (x : —r) 

b] 


dy 


the factor is « and the equation becomes 


8 
da Tad XG") 2 Dh ee ~ r) Q. 
dx 
whence 
ew’ U «2 oH. ~féel(x 2p 

z ee (xg iy °) rae eft ty ") 3 Le fark dy ?) 4(9), 
(y) being an arbitrary function of y. 

We shall take as an example the equation 


dz ‘ dz , 
x y— = Nz. 
de 7 dy 
dy 
Let — = dt, and therefore y = ¢«'. Then 
y 
dz 1 dz z 
ats a "2 
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The integral of which, by the preceding formula, is 


(ss @ 
a Sel it~ 3) (2), 


—log x s 
or z= &" log x dt ? (t) 
= a2" ¢o(t — log 2) = 2" y(log y — log 2) = ary (! z), 
; dz dz a2 
The equation ea~¥ rs - 


may be integrated in the same way. Or if we change both the 


_ 
independent variables, making i = dt and = = du, it becomes 


(= - aad 
du dt z 


d 


which gives z = oA fdwe " wit g2u --t + . * H9(0) 


tt 


i du eM eta i 4 (1) 


a =t 


—g- + ot + u) 


H 


24 
= + vay), 


as wu = log x, t = log y. 
dz 
The equation y z: 


C4 dz 
x 


+2?— =2 
dy 
may, by changing y dy into 4 d.y®, be transformed into 
dz d 1 
~ + (20 - —;)e=0 
The integral of which is 


ce? l 
DO a on day’ — -:) (y*) 
_ St 


=e 7 2 i iH 4(y2) 
—z? é Pmt f= 
=e ov ge dv! V) gy?) 
é oe 
= ode 3S TH g(yry 
d 
- , ay Selog (atv xttyt) | (yy?) 
d 
=e 4 Sa 4 V2? + y) o(y 
(x + y) o(y? — 2°). 
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Other equations, which at first sight do not appear to come under 
this form, may be reduced to it by a proper assumption of a new 
independent variable. For instance, the equation 


dz dz 
(@+y)7 + (y—2) =* 
is converted into 
dz Zz 
” de du 
by supposing y =u — 2. These however are particular cases 
which come under no general rule. 


= 2, 


D. F. G. 


X.—ON STABLE AND UNSTABLE EQUILIBRIUM.* 


Bae BP cesses be any number of forces acting upon a material 
system, and let dp, dp’, dp"...... be the virtual velocities of their 
points of application, estimated positive in the directions in which 
the forces act; then, when P, P’, P”...... are in equilibrium, the 
equilibrium is stable when {(Pdp + P’dp' + ...... ) is a maximum ; 
and unstable when this is not a maximum: and vice versd. 

We shall demonstrate this in the case of one rigid body; and the 
same reasoning can easily be applied when the system contains 
several rigid bodies. 

Whee PF, P, PO ...02. act upon a rigid body, they can always be 
reduced to two forces, but not to one. Let R and R’ be these 
forces, and dr, dr’ the virtual velocities of their points of applica- 
tion: then 


Pdp + P’dp' + ...... = Rdr + Rdr’, 
“.u = f(Pdp + P'dp’ ere ) = J (Radr oa R'dr’). 


Hence du = Radr + R'dr’, 
and d’u = Rd?@r + dRdr + R'd?r’ 4 dR‘dr’. 
Now when P, P’, P’...... are in equilibrium, 
Pdp + P’dp’ + P"dp" + ...... = 0, 


“. Rdr + R'dr = 0, 2. du = 0; 
and, since dr and dr’ are independent of each other, 
R = 0 and R' = 0, 
2. @u = dR.dr 4 dR'.dr’, 
When « is a maximum, du is negative, and therefore (since dr 
and dr’ are independent of each other) dR.dr and dR’.dr’ ave both 


* From a Correspondent. 
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negative. Hence if dr be positive, dR is negative; and vice versa: 
and so with dr’ and dR’. 

From this we learn, that if the body receive any slight displace- 
ment from the situation of equilibrium, two small forces dR and dR’ 
are brought into play, which act opposite to the directions in which 
the points, at which they act, have moved in consequence of the dis- 
placement of the body. Hence the effect of these forces is to draw 
the body back towards its position of equilibrium: and therefore the 
equilibrium was stable. 

Conversely, in order that the equilibrium may be stable, the forces 
dR, dR’ put in play by the displacement must act so as to carry the 
body back to its situation of equilibrium: and therefore the virtual 
velocities dr and dr’ must have different signs to dR and dl’ respec- 
tively: hence d?u is negative, and therefore wv a maximum. 

If uw be a minimum we may shew in the same manner that the 
equilibrium is unstable, and vice versa. 

If u be neither a maximum nor a minimum, then d2u (or d°u, if 
du vanishes) will be positive for some displacements, and negative 
for others; and will therefore be apparently dubious, though in fact 
it will be unstable; for although a displacement may be found which 
shall cause the body to return to its position of rest, yet in oscillat- 
ing through the position of rest it may come to one, from which it 
will not return to its position of rest. The converse of this is also 
true. 

Hence, if by stable equilibrium we mean, that for all small dis- 
placements the body oscillates about its position of rest, then the 
equilibrium is stable when wu is a maximum, and in no other case: 
in all other cases the equilibrium is wastable. 

N.B. We have supposed that Pdp + P’dp'+... is a perfect dif- 
ferential. If this be not the case, the above reasoning is still true ; 
but the result should be stated thus: when the equilibrium is stable, 
the first differential of Pdp + P’dp'+... (or if that and the second 
vanish, the third) is negative, and vice versa: in all other cases the 
equilibrium is unstable. 

If the system consist of several bodies, and R, R’, R’, R’’...... 
be the smallest number of forces to which the system of forces act- 
ing on the bodies can be reduced, then 

du = Rdr + R' dr’ + Rar" }- Re da” + 
in which dr, dr’, dr’, dr’” ...... are independent of each other: and 
the reasoning upon these and the forces R, R’, RY’, R"...... will be 
as above. 

The following is an application of this principle to fluid equi- 
librium. 


incompressible fluid, every particle of which is attracted towards 
a centre of force, varying as a function of the distance; and to 
determine whether the equilibrium be stable or unstable. 


Pros_tem. Required the form of equilibrium of a mass of 
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Let the fluid be referred to polar coordinates, the centre of 
force being the origin, @ the angle which 7 makes with a fixed line 
through the origin, and w the angle which the plane through r 
and the fixed line makes with a fixed plane passing through the 
fixed line: and let (as usual) cos 0 = p. It is evident that the 
fluid must completely surround the centre of force: let r and r’ 
be the radii of the external and internal surfaces; these are 
functions of » and w. The volume of the fluid is given: let 
ce be the radius of the sphere which has the same volume: hence 

4, 

: 


1 27 pt io a 29 E J 
; 38 = ae } 7 du dw dr = a/_ J, (r? — r) du dw. 


Since r* and r? are functions of » and w, they may be expanded in 
series of Laplace’s coefficients: let 
r= a + a (7, ns ¥, re )= a + ay, 
and of? — 2 4 0(¥, + Y, + ...... ) = a3 + ay’, 
where a is a numerical coefficient, the use of which will be seen 
hereafter. 
Then, by a fundamental property of these functions, 


al .27 
| an y du dw — 0; 
and similarly of y’. Hence we have, by the above equation, 
A&8— a — a’, 

Let f(r) be the force of attraction on a unit of mass at a 
distance r from the centre of force: then f(7) cm is the attraction 
on the element of the mass ém: this attraction acts in the line 7, and 
tends to shorten 7: hence, in this ease, 


v= [3 — f(r) dm? dr; and this s./$—f(r) dr? om. 

The symbol of integration = refers to the differential of the 
mass ém. Now ¢m equals (as above, if density = 1) 

r? sin 0 dd dw dr = — r* du dw dr, 
and the limits of 7, w, p are respectively r and r, O and 2z, 
1 and — 1: then, putting 
[—S(r) dr =f, (r), and fr f(r) dr = ¢ (7°), 

we have 


2 ¢ 1 27 os 
u =f" JS, i rf (r)dudwdr = /_ ¥ 3 $@ (13) — (r')} dp dw. 


But (1?) — @ (1%) = o(a® + ay) — g(a? + ay’) 
= $(a°) — 9(a2) + ay ¢ (a) — ay 9 (a) 

+ $a°y? og’ (a) — da’y’? o" (a) +... 
in which the accents to g denote differentiation with respect to 
a® and a*. Hence 

u = 4r@ (a?) — 479 (a’*) 


+ 4a?/ m3 "$y2o (a®) — y'29" (a du dw 4 
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‘ ; / . » . 
It is very easily shewn, that »’ (a*) = — 5 7 J(@), and is there- 
Ya* 
fore negative in every case of attraction, (but positive if the force 
were repulsive). 

From these calculations we gather the following results. If a be 
a very small quantity, the increment of u above what it becomes 
when a = 0 (that is, when r =a and r = a’) does not involve 
the first power of a: hence, when r =a and r' =a’, we have 
du = 0, and therefore the fluid is in equilibrium. 

Hence a form of equilibrium is that of a hollow sphere of any 
dimensions (so long as the volume be constant), and having the 
centres of the bounding surfaces at the centre of force. 

And the equilibrium of the external surface is stable, because 
for all values of y when y'=0 the multiplier of a® is negative, and 
therefore uw is a maximum: and the equilibrium of the internal 
surface is unstable, because for all values of y' when y = 0 the 
multiplier of a* is positive, and therefore u is a minimum. On 
the whole, the equilibrium is dubious. 


Cor. 1. If a’=0, then y’=0, and the equilibrium is altogether 
stable. 

Cor. 2. If the force were repulsive, and the fluid were held 
in equilibrium by a rigid spherical envelope; then y = 0, and the 
multiplier of a* would be negative, since @”(a*) is then positive, 
and u is a maximum, and the equilibrium of the internal surface is 
stable. 

J. tA, Pe 


XI.—ANALYTICAL GEOMETRY OF THREE DIMENSIONS. 
No. II. 


WE proceed to give some instances of the advantage of symmetry 
in applying the Differential Calculus to Geometry of three dimen- 
sions. The following form of the equation to the tangent plane is 
probably known to many of our readers; but, as Leroy has only 
deduced it from that in terms of the partial differential coefficients 
of one of the coordinates with respect to the other two, we will give 
an independent and easy proof of it. 


1. If F(a, y, z) = 0, be the equation to a surface, the locus of 
the tangent lines drawn to it at a point (a, y, z) is a plane whose 
equation is 
_— dF ett «<2 d¥ ; dF 

"” dx ' we dy 7 
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The equation to any line drawn through the point 2, y, 2, is 


‘ ‘ / 
Oe oe ne icici 





l m n 


This line will in general be cut by the surface in one or more 
other points. Let the coordinates of the nearest of these be 2 Yy» Zp 
and let each member of equation ( (2), when 2’, ’, 2’ become 2) Yp Zp 
be assumed equal to 7, then 


a=ath, y=ytmy, 42 +07; 
but F(a, ¥ s,) = 0; 
therefore F $x 4+ ly, y + mr, z+ nr} = 0. 
Expanding and observing that F(a, y, z) = 0, 
dF : F d F d2F 
i no = BH BCe ) Ae ceceee == 
.* _-* dz +5(S + c.) + 


When the es becomes a tangent, the points (a, y, 2), (@p Yv 2») 
approach indefinitely near to one another, and r becomes indefinitely 
small ; therefore, taking the limit of the preceding equation, 

dF iF iF 
i +-m— +2—=0, 
da "dy dz 

Multiplying the several terms of this by the several members of 
equation (2), we obtain 

dF 
+ (2 — 2) = == 0, 


for the equation to the locus of the tangent lines. 


d¥ 
@ - 2) + (y'- ”) ty 


2. Since the normal line is perpendicular to the tangent plane, 
its equations are 
a—e y-y zz 
en a = gp 
dx dy dz 


3. We shall next investigate an expression in terms of the partial 
differential coefficients of F (2, y, z), for the radius of curvature of 
a section of a surface made by any plane passing through the nor- 
mal at a given point. 

Let dF =U, a¥F = V, aI = W, then the equations to the 

dz dy dz 
normal are 








ws _¥-9 2-8 

— Vv Woo (4). 

Again, if ds be an element of the section, and dz, dy, dz its 
dz 


dy , oe 
, oe i are the cosines of the angles which it 


projections, 
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makes with the axes, and therefore the equation to the normal 
plane to the section is 
Ce av) dx + (y/ — y)dy + (2 — 2)de = 0.0855 (5). 
The line of intersection of two consecutive normal planes will be 
determined by (5) and its differential, which is 
(x —2x)d*x + (y' —y) @y + (2x —2z) dz —ds? = 0...... (6); 
and the intersection of this line with the normal to the surface will 
evidently be the centre of curvature of the section. We may sup- 
pose 2’, y’, 2’ to belong to this point; they will then be the same in 
(4) as in (5) and (6): and assuming each member of (4) equal to 
Q, and substituting for 2’, y’, z' in (6), we have 
(Udret+Vdy+W d’*z) Q — ds* = O. 
But if p be the radius of curvature 
pe =(e# — 2) +(y — yy) + (2 — zy =(U* + Vt + WW?) Q?, 
therefore 
VU2 + V2 + W?2.ds? (7) 
U dx } V dy + W dz eee eee eee i . 
This expression may be transformed as follows. Since 


Udx + Vdy + Wdz = 0, 


p= 


ee ae , d°F d’F d°F 
differentiating, and assuming ——~ = u, — = vy, — = w, 
dx dy? dz* 
d°*F _ atF - dF ; 
= wu, =v, = -, 
dy dz dz dx dx dy 
U dx 4+ Vdty + Wd?z 4+ uda® 4 vdy®? + wdz’ 
+ 2u' dy dz + 2v dz dx + 2w' dx dy = 0. 
é : hi ti . dx dy dz 
‘ ying this equ: , and assuming =/,+-= — =f, 
Employing this equation, and assuming = z= = 


(7) is changed into 
7 VU2 4 V2 4 W? (8) 
p= Pu + my + n2w + 2mnu 4 Qnle’ + 2lmw' , 
4. It may be as well here to show how formule involving the 
partial differential coefficients of F, may readily be transformed into 
others in terms of the partial differential coefficients of one of the 
coordinates as z. Suppose the equation to be put under the form 
f(a, y) - z = 0, so that F(a, y, z) = f(a, y) —2. Then 
dF d dz dF dz dF ) 
=" = a f(a, y) “_ b ] hi oR = - 2 , “Omg = —_— i. 
dr da* ‘ dx dy dy’ dz 
°F 22 °F Itz °F 
5 Al A ll = 0, = 
dx’ dx dy? dy® dz* 
1°F d°F d°F d*z 
+l = 0, —— 0, —<<—< = ——, 
dy dz dz da dx dy dx dy’ J 
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If we substitute these values in (8) and adopt the usual notation, 
we obtain the known expression 


meee ns f 7 4 
P= Re + 2lms + mt’ 


5. To find the sections whose radius of curvature is a maximum 
or minimum, and the values of those radii. 


In formula (8) let += YU* + V* + W* = P, then 
: = Pu + mv + n?w + 2mnu' + Qalr’ + Zlmw’...... (10), 
the quantities /, m, m being connected by the two equations 
Bide a? > of oe D nccccsscvech ID) 
UW + mV + aW = 0............(12). 
Differentiating these equations, 
(lu + mw’ + nv’) dl + (le 4+ mv + nu’) dm 
+ (lo' + miu! + nw) dn = 0......... (13), 
ldi + mdm + ndn = 0.......:. (14), 
Udl 4+ Vdm 4+ Wdn = 0........ (15), 
(13) + (14) x A + (15) p gives, on equating to zero, the coefii- 
cients of each differential, 
le 4+ mw’ + nv + Al + pU = 0 
lo’ + mv + ni’ + dm + pV = 0 | re (16). 
le’ + mu 4+ nw + An 4 pW 0 


Multiplying the first of these by I, the second by m, and the 
third by n, adding, and reducing by C10), CLT), (12), 


P 
+rA =O 
p 
Substituting for \ in (16), they become 
fe P\ / U Y 7 
(w — J f+ wm 4 vn + pU =0 
P 
4 ; 
wl 4 (c — = | m -| un + pV == 0}...(17). 
p 
P 
vl + “um + (w — -) n+ wW =0 
Pp 





To obtain the equation for p, it remains to eliminate /, m, x, Ht 
from equations (17) and (12). The result is 
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>) » >» 
U2 (« _ “) (w -- ~) + V? (w - -) (u - ~) 
p p p p 
> ») 
— 2VWu (u *) — 2WUv’ (c *) — 2UVw' (w _ -) 
p p p 


— U2u'2 — V2v'2 — W2w'2 
4+ 2VWe'w' + 2WUw'u' + 2UVu'e' 
6. If the equation to the surface be of the form 
tH) (x) + x (y) oa Ww (z) = 0, 
which includes, among others, all the surfaces of the second order 
when referred to their principal diametral planes; U, V, W are 
functions of 2, y, z alone, respectively, so that wu’, v', w’ are all 
zero, and the equation (18) simplifies considerably. In short, in 
this case, equations (17) give immediately 





1U 1V iW 
i=—,, m = ———,, n= ,; 
“u—— v—-~ w—— 
p p p 
and by substituting these values in (12), we have 
U? _ Ww? 
p+ p+ p WG cissvisicch 
“u—— v-- w-- 
p p p 


7. Exampce. To find the principal radii of curvature of an 
ellipsoid. 


In this case 


: x y” 2° 
F(a,y%z)= [+E +5 -1=0, 
a’ be ce 


therefore 
Or o° O- (9) a) a 
a. , , @ 7) 2 2 
U = ——» \ =F? W =, u =-s» o= = w=, 
a’ b c a® b? c 


and the value of P is 
x yy" 2 2 
9 _ A _— 
a al Ba ae 
c pP 
if p denote the perpendicular from the centre on the tangent plane. 
Substituting these values in (19), we have 


9 


x? y 


a (pp — a") * 8 (pp — 0) * & (gp — ct) ~ 
which may easily be put in the form 
24202 
pr — far 4 4 c&@— (22 + 4 eel 4 % ts 


P p* 











of 
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8. As another example, we propose to find the principal radii 
of curvature at any point of the surface whose equation is 
ata? 4 by? 4 cz? = (x? 4 y? 4 2%)? ...... (a), 
and which is the locus of the extremity of the perpendicular from 
the centre on the tangent plane to an ellipsoid. Put 
a? 4 y2 22 = 72, 
then F (a, y, z) = a*x? + By? + cz? — 4, 

U = 2x(a? — 2r*), V = 2y(b?— 2r?), W = 2z(c? — 2r*), 
u=2(a?—2r?—42"), v= 2 (b? —2r?—4y"), w= 2 (c?~ 2r?— 422), 
ul’ = — 8yz, v' = — 82a, w' = — Bry. 

Instead of substituting these values in equation (18), we shall 
employ the equations from which that is deduced, namely (12) and 
(17), because by this method it is more easy to obtain a result of 
a simple form. Of these, (12) becomes 

Ix (a? — 2r?) + my (0? — 2r?) + nz (ce? — 2r?) = 0... (b), 
and the first of (17), 


P 
l (w — 27? — 4a? — “) — 4maxy — 4nxz + px (a? — 2r?) = 0, 
P 


) 


I 
or l (« — 2r? — -) — 4a (lx + my + nz) + pa (a? — 27?) = 0. 
p 
Assume be - my 4 nz = Re .....cseeseeee (€)3 


then we have 


Pp 
1 (a — 2r2 — - ) — 4hxe + pu (a? — 2r?) = 0} 
Pp 


) 
m (s — Qr? — -) — thy + py (8? — 2r?) = Ob... (ad), 





P 
n ( — 2 — =) — 4hz + pz (c? — 2?) = 0 
Pp J 
and we shall obtain the equation for p by eliminating /, m, n, hk, My 
from the five equations (6), (c), (d). For this purpose, multiply 
equations (d) by 2, y, z respectively, and add them, reducing by 
(a), (6), (ce), then we find 
P 
k ( + 1) + pr? = 0. 
Pp 
By this relation, eliminate » from equations (d), and we have 
_P 
4 (a2 — 72) 72 2 972) — 
(a r?) 72 + (a r?) al 
P "gf? 
Pp 


with corresponding expressions for m and x. Substituting them in 
(6), and denoting 


F (a, x) + F(b, y) + F(e,z) by =F (a, x), 


‘= 





a* — 2r? — 
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» 
4 (a? — 7?) (a2 — Qr*) r? + (a? — ony 





y p C . 2 = 0; 
a? — Qr? — — 
pP 
but by (a) x(a? — 7) a? = 0: 


multiplying the last equation by 47°, subtracting from the preceding, 
J 
and dividing the result by — , 
p 
™ a‘z* 
a ——— == 0, 


2 : 
which is the simplest forni in which the equation can be presented. 
The value of P will be found to be 
2 Vaix? + by? 4 e822, 
The length of the perpendicular from the centre on the tangent 
plane to this surface is 
ri 
V ax? 4 by? 4 e422’ 
Ort 
therefore if p denote that line, P - _ . 
Pp 
9. The directions of the sections of greatest and least curvature 
will be found by eliminating \ and p» from equations (16). The 
result is 
(lu + mw! 4 nv’) (mW — nV) + (lw + mv + nu’) (nU — IW) 
+ (lv + mu’ 4- nw’) (IV — mU) = 0, 
or (Ve — Ww') ? 4+ (Ww! — Un’) m? 4 (Ua — Vo’) n? 
as sWe' —Vw'+U (v—w)} mn + s$Uw' — Wu! + V (w—u)} nl 
t $Vu'—Ur' + W (u—v)} lm = 0 ...... (20), 
which equation, together with (11) and (12), determines the values 
of 1, m,n. It is easy to see that there will be two sets of values of 
2. m*, n®, and no more. 
10. To prove that the directions of greatest and least curvature 
are at right angles. 
For abbreviation, write the equation (20) 
LP 4+ Mm? + Nn? + L'mn + M’nl + Nim = 0. 
Substituting the value of 2 derived from equation (12), 
(LE + Mm* + Nim) W? — (L'm + M7) (JU + mV) W 
+ NU 4+ mV)* = 0, 
or (MW? + NV? — L'VW) m? + S2NUV— W(L'U + M'V)} lm 
+ (NU24LW2-MWU)2—0, 
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: ° — . f 
a quadratic equation in —; whence we find, that if /,, 2, are the 


l 
two values of /, and m,, m, those of m, 
mm, NU? + LW? — M'WU 
ii, MW? + NV? — L’'vw’ 
so that if we assume 
Ll, = K(MW* + NV* — L'VW), 
we shall have mm, = K (NU* + LW* — M’WU), 
and also nn, = K(LV* + MU? - N'UV). 
Hence, the cosine of the angle between the two directions, or 
Ll, 4+ mm, + nn, = K 3L (V? + W?) + M(W? 4 U?) 
+ N (U2 + V?) — L'VW — M'WU — N'UV3. 
The value of the quantity within the brackets is zero, as may 
easily be seen by recurring to equation (20) for the values of 
LM, &e, and collecting the terms multiplied by w’, 2, &c. respec- 
tively. Therefore the two directions are at right angles. 
11. To find the lines of curvature of a given surface, or the 
direction in which we must pass from one point to a consecutive, in 
order that the normals at the two points may meet. 


Retaining the same notation, the equations to the normal at 
a point a, y, 2 are 
aaa y-—y_2-? 
U ~ Vv Ww 
Let each member be assumed equal to Q, then 
a = 2 + QU, y¥ =y + QV, 2 =z + QW. 

If 2’, y', z belong to the point in which two consecutive normals 
meet, the differentials of these three equations must be verified at 
the same time with them, therefore 

dx + QdU + UdQ = 0, 
dy + QdV + VdQ = 0, -............ (21), 
dz -f QdW | WdQ =< ©. 


or udx + wdy + vdz + : dx + U aa =O; 
Q Q 

wdx + vdy + udz + G dy+ V z = 0, 

vdx + udy + wdz + a dz + W * = 0. 


These equations are of the same form as equations (16), therefore 
the elimination of Q and dQ from these must give the same relation 
between dx, dy, dz, as the elimination of \ and p» from (16) gave 
between /, m, n, namely (20). And if 
da® + dy’ + dz? = ds’, 
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ll dy e will be the cosines of the angles which an element 
ds ds’ ds 

of a line of curvature makes with the axes, and they satisfy the 
same system of equations as J, m, n, namely (11), (12), (20). 
Therefore the directions in which two consecutive normals meet 
coincide with those of greatest and least curvature. 

For the sake of shortness, we shall eliminate Q and dQ from 
equations (21), instead of their expanded forms. This is easily 
done by cross-multiplication, and the result is 

U (dVdz — dWdy) + V (dWdx — dUdz) 

+ W (dUdy — dVdx) = 0......... (22): 
from this equation and that to the surface, the differential to the 
projection of a line of curvature on any coordinate plane may be 
found. 

12. Exampce. To find the lines of curvature of an ellipsoid. 
Let the equation be 


Fi y 2 
—4+74+5=1; 
a b c 
Qu : 2 , 22 
then U = , v= 4 W=-—, 
a’ b° c 
2 2 4 
dU =—dz, dV=, dy, dW = —dz. 
a” Bb ce 


Substituting in (22), 
(b? — c*) w dy dz + (ce? — a) y dz dx + (a? — *) zdxdy = 0. 


~ 
~ 


Multiply by =, and substitute the values 
r Cc 


4 a yy 2zdz x dx y dy 
prb-o-} wae 
therefore 
d 
10 — Ay ady + (2 — yan (2 4 24) 
@ 


x2 7 
— (a2 — ) (1 ace? ae y) de dy = 0, 


S| 


or, reducing, 





b? — c? 
ay dy’ 
@—c | ec . ae at ct 
+ f —— 2° — —____ y* — (a* — 5° | dx dy — —— ay dx* = 0. 
| a b° ; J ‘ a? ‘ 
To integrate this, multiply by 4ay, and let 2* = u, y® =v, then 
6b? — ec z 
= u dv? 
a*—c? 82 — @? ; a— 
t — eae , v—(e& By | du dv — : a du? = O, 
| a b | a? 
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or (“E dv + a = du) (u dv—v du) —(a°—b*) du dv=0, 
a 





b’ 
as du 
(a — b*) - 
du dv 
orw—v—=F . . 
dv b2 — ce a? — c2 du’ 
b a dv 


an equation of Clairaut’s form, whose solution is 


C (a? — 3?) 








oe se 6? — c? , a? — ¢?’ 
b +s a? 
oe ae C (a? - b) 
b — c? ,a@ — cr’ 
Bb iis a 


the equation to an ellipse or hyperbola. For the method of deter- 
mining the constant, see Leroy, No. 422. 


XI1.—MATHEMATICAL NOTES. 


1. Taylor's Theorem. The following proof of Taylor’s Theorem 
by means of the separation of symbols of operation from those of 
quantity, may be interesting to some of our readers. Let there be 
an operation, which we may for convenience call D, which per- 
formed h times on f(a) converts it into f(a +h); that is, let 
D' f(x) = f(x +h), we wish to determine the nature of D. 
fet 5-9) 

h 


: d er 
Now we know that s f(~) = when fh is in- 
dx* 


definitely small, therefore 
@ 2 D' f(v7) — f(x) _ D'— 1 
ded (z) = h a h F(®) 
when f is indefinitely small. Therefore, taking the symbols of 
operation separate, 
d PDr—1 


toda when 4 — 0, 
dx t 


1 
d\h 
which gives D = (1 +h x) when h = 0. 
dx. 
Expanding by the Binomial Theorem, and making 4 = 0, we find 
1? 1 d8 
D=—1 -f- @ b 4 : 


dx 1.2 dx 1.2.3 dx 


+ &e. 
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qd 
or D — ¢47; 
ae 
therefore f(x + h) = « 4 f(x), 


which is Taylor’s Theorem. 


2. Problem from the Papers of 1835. If p, r be the perpen- 
dicular from the origin on the tangent plane and the radius vector 


. J p> . . 
of any surface, then Pp will be the perpendicular on the tangent 
> 


plane at the corresponding point of the surface, which is the locus 
of the extremity of p. 

Let 2, y, z be the coordinates of the first surface, a, 6, y of the 
second, As p is the perpendicular on the tangent plane, the equa- 
tion to the plane is 

ax + By + y2 Pp a 4 Pp? 4 y? 

Since this plane is a tangent to the surface, this equation will 
hold good if for x, y, z we put a+ dz, y + dy, z +dz. Whence 
adv 4+ pdy + Y ae =D (2). 

Now if V =0 be the equation to the surface of which a, }, y 
are the coordinates, and P be the perpendicular on the tangent 
plane ; then 
dV + , av ‘ dV 
da "dp * dy 


(av 2 dV\? dV a’ 
re ye 
\ da 3 dy] | 
1V 1V 1V 
also, —— da + aa dp + : 
da dj3 dy 
Now differentiating (1), considering a, 3, y, x, y, z as variables, 
and paying regard to (2), we have 
(w — 2a) da 4+ (y — 28) dB 4 (2 — 2y) dy = 0 
A (3)-(4) gives, on equating to 0, the coefficients of each of the 
differentials 
dV dV dV 
>» ss 


— xa— 2a, Xr = 9 25, OX —2z— 2. 
da dp . . 


a 





ee 


: dy 
Substituting these in the expression for P, it becomes 
p— __ 2 (@? + + 7°) — (aw + fy + 72) 
py ees tla t P+ — (att py + ARE 
which by (1) is reduced to 
“et a? + 6? ar Pr 


t 





2 








